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The quantum analog of the joint probability distributions describing a classical stochastic process 
is introduced. A prescription is given for constructing the quantum distribution associated with a 
sequence of measurements. For the case of quantum Brownian motion this prescription is illustrated 
with a number of explicit examples. In particular it is shown how the prescription can be extended 
in the form of a general formula for the Wigner function of a Brownian particle entangled with a 
heat bath.. 

I. INTRODUCTION 

The notion of joint probability distribution is basic to the description of classical stochastic processes. The purpose 
here is to describe the extension to the quantum regime, giving the prescription for constructing the quantum joint 
probability distribution associated with a sequence of measurements. The prescription is illustrated with a number 
of examples. These are an important part of this work, since they show how the prescription can be used to calculate 
a variety of quantities of practical interest. 

The idea is simple: The system is initially in thermal equilibrium. A first measurement prepares a state that then 
develops in time according to the underlying dynamics. Then a second measurement prepares a new state. And so 
on until the final measurement in the sequence. 

In this connection it is necessary to consider the description of quantum measurement. In its most naive form, 
found in many textbooks, a quantum measurement of a dynamical variable is described as a projection of the system 
into an eigenstate of the variable, with no memory of the previous state. In particular for a variable with a continuous 
spectrum, such as the position of a Brownian particle, with no square integrable eigenstate, this naive description 
is unsatisfactory. In an earlier publication jointly with J. T. Lewis the description of measurement as applied 
to quantum stochastic processes was addressed in some detail. Since that earlier publication may not be accessible 
to all readers, in Section [ni] a summary is given of the essential features of quantum measurement as they apply to 
the definition of the distribution functions. The reader will observe that the description given there involves no new 
theory of quantum measurement. Rather, a prescription is adopted based on that used by many authors making 
practical calculations related to real experiments. 

The plan of the paper is as follows. To begin, in a brief Section |TT] the joint distribution functions of classical 
mechanics are described. The quantum joint distribution functions, which are a close analog of the classical quantities, 
are introduced in Section IIIII . There the key result is the prescription p.29p for the joint distribution function 
associated with n successive measurements. For the case of quantum Brownian motion this prescription can be 
readily evaluated to give explicit closed form expressions. Therefore in Appendix [Xl a review is given of those aspects 
of the theory of quantum Brownian motion that will be useful in the applications. There the key quantities needed 
for the later discussion are the commutator and the mean square displacement, given by the general expressions (|A7p 
and (|A11[) . Later in Appendix [X] these expressions are evaluated explicitly for the Ohmic and single relaxation time 
models, and the results compared with approximate expressions obtained by master equation methods. In Section HVl 
the results of Appendix [X] are used to evaluate the characteristic function associated with the distribution function 
describing n successive measurements. There the important result is the expression (j4.2p for the characteristic function, 
where it is seen explicitly how in classical mechanics, where the commutator vanishes, the effects of measurement 
can be separated from the underlying stochastic process. As an application of this result, an explicit expression is 
constructed for the pair distribution function associated with wave packet spreading. In Section |V] we discuss the 
probability distribution, which corresponds to what in elementary quantum mechanics is "the square of the wave 
function" . A general expression is obtained that is illustrated first with the example of wave packet spreading and 
then with the example of a "Schr odinger cat" state. In either case the discussion includes the case of a free particle 
as well as that of a particle in a harmonic well. Finally, in Section IVII the Wigner function is introduced. There the 
key result is the simple formula (|6.5p for the Wigner characteristic function (the Fourier transform of the Wigner 
function). The Wigner function corresponds to the phase space distribution of classical mechanics, but is definitely 
not a probability distribution (it is seen explicitly in the examples that the Wigner function need not be positive) 
and cannot be the result of direct quantum measurement. Nevertheless, the Wigner function is useful for describing 
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the results of measurement. In particular the probability distribution in coordinate or momentum are obtained by 
integration over the conjugate variable. 



II. DISTRIBUTION FUNCTIONS IN CLASSICAL BROWNIAN MOTION 



In the theory of classical Brownian motion a stochastic process is completely described by a hierarchy of probability 
distributions. Here the standard reference for the physicist is the review article by Wang and Uhlenbeck reprinted 
in the "Noisebook" [3]. For a stochastic variable y{t) one introduces 

W{yi,ti)dyi — probability of finding y{ti) in the interval dyi about yi, 
W{yi,ti\y2,t2)dyidy2 = probability of finding y{ti) in the interval dyi about yi 

and y{t2) in the interval dy2 about y2, 

and so on. (2-1) 
This hierarchy must satisfy the following more or less obvious conditions 

1. Positivity 

W{yi,ti;y2,t2; ■ ■ ■ > 0. 

2. Symmetry 

W{yi,ti] 2/2, ^2; • • ■ ; y-n, tn) is a symmetric function of the set of variables yi, ^i; y2, ^2, • • • i Vn, tn- 

3. Consistency 

W{yi,ti;y2,t2; ■ ■ ■ = J dyn+iW{yi,ti;y2,t2; ■ ■ ■ ; y„, t„; 2/„+i, t„+i). 

Note that consistency corresponds to conservation of total probability, 



1 = 



J dyi J dy2--- J dynW{yi,ti;y2,t2; ■ ■ ■ (2.2) 



A theorem of Kolomogorov states that if the hierarchy satisfies these conditions there must exist an underlying classical 
process Q. That is, there must exist an ensemble of time-tracks y{t) such that the W^s are the weighted fraction of 
time tracks that go through the appropriate intervals. 

A natural question is How is this description changed in the quantum case? The answer will be seen in the following 
Section, but for now one can say that the essential change is that in the quantum case the symmetry condition no 
longer holds. 



III. QUANTUM DISTRIBUTION FUNCTIONS 



The system considered is that of a Brownian particle coupled to a heat bath, a system with an infinite number of 
degrees of freedom. The quantum mechanical motion of this system is described by a microscopic Hamiltonian H 
and corresponds to a unitary transformation of states in Hilbert space, 

- C/(i)*(0), (3.1) 

where ^(i) is the state vector at time t and 

U{t) ^ exp{-iHt/h}. (3.2) 

However, one does not have precise knowledge of the initial state. Instead, there is an initial density matrix. The 
density matrix is defined as an operator p{t) in Hilbert space such that (<I',p(i)$) / ($,<&) is the relative probability 
at time t that the system is in any given state $ [5|, 16j . Note that consistent with this definition p must be a positive 
definite Hermitian operator. Its time development follows from (j3.ip . 

p{t) = U{t)p{0)uHt). (3.3) 

Before introducing the distribution functions, it is necessary to make some general remarks about measurement in 
quantum mechanics. By "measurement" here is meant "measurement with selection" (what Pauli in his famous Ziffer 
9 called "measurement of the second kind" Q) so that measurement irreversibly changes the state of the system. 
When discussed in general terms in textbooks, the accepted description of this change of state is framed in terms of 
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measurement of a discrete variable (Hermitian operator with a pure point spectrum). Let B be such a variable, with 
b an eigenvalue and its associated projection operator, so that 



B = Y^ bPb. (3.4) 



Then the effect of a measurement at time ti whose result is that the eigenvalue b is in the interval M is to instanta- 
neously transform the density matrix, 

Pih) ^ PMp{tl)PM, (3.5) 

where Pm is the projection operator associated with the interval, 

Pm = (3.6) 

be A/ 

Here instantaneous means that the duration of the measurement is short compared with the natural periods of 
the system. The prescription p. 51 ) may be obtained from various assumptions about the optimal character of the 
measurement (such that the disturbance of the state is somehow minimal). See, e.g., Liiders Goldberger and 
Watson [9|] , Furry , or Davies and Lewis [ll| . 

But the prescription (|3.5p is too restricted for our purpose; it represents too limited a class of measurements. We 
must consider measurements of limited precision and involving operators with a continuous spectrum, such as the 
position of a Brownian particle. For guidance as to how to generalize the prescription, we look to such practical 
fields as the theory of angular correlations (see, e. g., the article by Frauenfelder and StefFan §3) or the theory of 
polarization in multiple scattering (see, e. g., Wolfenstein ^1§\ §4). There one associates a transition operator T (in the 
scattering case this would be the Wigner T-matrix) with a measurement with a given result (e.g., the observation of 
an emitted gamma ray in a given direction) and represents the transformation of p brought about by the measurement 
at time ti by 

Pih) Tpit,)Tl (3.7) 

Note that this is the most general transformation that preserves the positivity of the density matrix. In the special 
case where the transition operator is a projection operator, one recovers the prescription (j3.5l ). To be consistent with 
the probabilistic interpretation of the density matrix, the transition operator must satisfy the general requirement 

||T||' = max^^Jl^<L (3.8) 

The diagonal matrix elements of p formed with respect to a complete set of stated are interpreted as the probabilities 
of finding the system in the corresponding states. The sum of these probabilities over all states is the trace and, as a 
consequence of the requirement p.8[) . this is reduced by measurement. Thus, 

Tr{rprt} ^ Tr{prtr} < Tt{p} \\Tf < Tt{p}. (3.9) 

In fact the ratio Tr{pT'^T} /Ti{p} can be interpreted as the probability that measurement will produce the given 
result. 

This reduction of the sum over all states of the probability of finding the system in each state is not a unique feature 
of quantum probability. In classical probability, where the probabilities after measurement would be interpreted as 
joint probabilities of the result of the measurement and of finding the system in the state, the sum of probabilities is 
also reduced. The difference is that in the classical case none of the individual probabilities will be increased, while 
in the quantum case some may increase. 

As a simple example illustrating all this, consider a spin 1/2 system initially polarized in the +z direction. This 
can be represented by the 2x2 density matrix 

Po=(ll)- (3.10) 



Note that the diagonal elements, which are, respectively, the probability that the spin is in the +z and —z directions, 
add up to 1. This is in accord with the convention that the trace of the density matrix is normalized to 1 prior to the 
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first measurement. Suppose that a measurement is made, for example by a Stern-Gerlach apparatus, whose result is 
that tlie spin is in the +x direction. Tlie transition matrix corresponding to this result is 

T-l{\\). (3.n, 

in this case a projection operator. The density matrix after the measurement is 

= TpoTt = Q I ^ . (3.12) 

The probability that the spin is in the +z direction has decreased from 1 to j while the probability that the spin 
is in the —z direction has increased from to j . Nevertheless the sum of the probabilities is ^, less than the sum 
before the measurement, One might ask: what happened to the probability, how is it that the probabilities after the 
measurement don't add up to 1? The answer is that there is another possible result of the measurement: the spin is in 
the —X direction. Repeating the above argument for this case, we find for the density matrix after the measurement 

= K ) • (3.13) 

Again the probability that the spin after the measurement is in the ±z direction is j. Thus, in either case the sum 
of the probabilities is the probability that the measurement produces the given result. These probabilities add up 
to 1, so overall probability is conserved. 

With these remarks as a guide, consider the measurement of a dynamical variable y. Assume that y is a variable, 
such as the position or velocity of the Brownian particle, with a continuous spectrum over all real values. In this case 
one can associate with the measurement a function f{yi) such that (here j/i is a c- number) 

(3.14) 

is the conditional probability that if the system is in state 4> the instrument will read in the interval dyi about yi . 
Here the choice that only the difference y — yi appears is made for convenience, since by choosing /(j/i) to be peaked 
the requirement that the measured value be somehow close to the actual value can be satisfied in an obvious way . If 
this conditional probability is to be normalized, one must require (/ need not be real) 

dyi\f{yi)f = 1. (3.15) 

An example to keep in mind is that of a "Gaussian instrument" [l^ . for which 



(2^2)1/4 em 4^2- 



where (Ti is the experimental width. The result of a measurement at time ti in which the instrument reads in the 
interval Ii is therefore to instantaneously transform p, 

pih)^ [ dy,f{y-yi)pih)f{y-yi)l (3.17) 



As remarked above, an instantaneous measurement is to be understood as one whose duration is short compared with 
the natural periods of the motion. 

The distribution functions are now constructed as follows. To begin assume that a.t t — (or in the distant past) 
the system is in equilibrium at temperature T, 

Pi^) = Po= rj.^^^_H/kTy (3-18) 
so p is initially normalized. If y is to be measured at a later time ti the system must move in time from to ti, 

Po^Uih)poU\h), (3.19) 
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according to (|3.3p . At ti a measurement is made, 

Uiti)poU\ti) ^ f dyifiy-yi)Uih)poU\h)fiy-yi)\ (3.20) 

according to (|3.17p . Finally, the system moves from ti to t, 

I dyJiy - yi)U{h)poU\h),f{y - y^^ 
Jh 

^ / dy,Uit-h)f{y^y,)Uih)poUHh)f{y-y,)^U^{t-t,). (3.21) 
Jh 

This last expression can be simplified somewhat by introducing the time-dependent variable (Heisenberg representa- 
tion), 

y{h) = UHti)yU{h). (3.22) 
Noting that U{t - h) = U{t)U{-ti) = U{t)W{ti), the final density matrix ((3?2T|) can be written 

/ dyiU{t)f[y{h)-y,]pof[y{h)~yi]^UHt). (3.23) 
Jh 

The probability that the measurement of y is in /lis the trace of this final density matrix. This same probability is 
interpreted as the integral over the one-point distribution W{yi,ti) over the interval Ji, so that 

/ dyiW{yuh)^ I dyiTY{U[t)f[y{h)~y^]pof[y{h)^yi\^U\t)}. (3.24) 
Jii Jh 

Since Ii is arbitrary, one can identify 

W{yi,h) = Tr{Uit)f[yih) - yi]pof[yih) - yi]^UHt)}. (3.25) 

Finally, the trace is invariant under cyclic permutation of the factors, so one can write 

Wiyi,h) = Tr{/[y(ii) - yi]pof[yih) ~ yi]^. (3.26) 

In the same way one can show that the two-point distribution is 

W{yi,h;y2, h) = Tr{/[y(i2) - y2]/[y(ii) - yi]pof[y{ti) - 2/i] ^[^(^2) - ^2]^} (3.27) 

and in general, using an obvious shorthand notation, 

W{1, ■■■ ,n) = Tr{/(n) • • • /(l)po/(l)t • • • f{n)^. (3.28) 

Finally, note that under cyclic permutation of the factors in the trace, one can write the expression for the rt-point 
distribution in the compact form: 

Wil, ■■■,n) = (/(l)t . . . /(n)t/(n) • • • /(I)) , (3.29) 

where the angular brackets indicate the thermal equilibrium expectation. That is, for a given operator O, 

{O) = Tt{Opo}. (3.30) 

The expression p. 2911 is the key result of this section. 

In connection with the expression ()3.29|) for the joint probability distribution, it should first be emphasized that 
it is time-ordered: < ii < t2 < ■ ■ ■ < tn- Also we should point out that in our shorthand notation the label 
applies to all the parameters of the measurement. Thus, for example the label "j" represents not only the value yj 
and the time tj but also the instrumental parameters such as the width aj . The symmetry property of the classical 
stochastic process refers to symmetry under permutations of the labels. In other words, the probability distributions 
of a classical stochastic process are symmetric under the interchange of all the parameters of the measurements. 
An inspection of the expression p.29[) for the quantum probability distribution shows that it does not have that 
symmetry, because the /'s at different times do not in general commute. However the quantum distributions still 
have the consistency property, providing the integral is over the results of the last measurement. This is sometimes 
called marginal consistency. 
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IV. CHARACTERISTIC FUNCTIONS 



In this section we consider the evaluation of the distribution functions when the dynamical variable y(t) is taken 
to be the position operator x{t) for quantum Brownian motion, introduced in Appendix \^ In evaluating these 
distribution functions, it is convenient in analogy with the classical case to introduce the corresponding characteristic 
functions, defined by 

/OO fOO ^ 

dxi--- / dxnW{l, ■ ■ ■ , n) exp{z fcjXj}. (4.1) 
-OO J —OQ j — 1 

An important result is that for quantum Brownian motion the characteristic function can be cast in the form 

air-- .n)=Kil,--- ,n)(^exp{iJ2kjx{t,)}^ , (4.2) 

where the factor ^(1, • • • , n) is given by 

Kil,...,n) = f[rdx,nx,^ ± ± (4.3) 

j=i''-°° 1=3+1 

In this expression the sums are to be taken to be zero when j = n. The importance of this result lies first of all in 
the fact that the effects of measurement are completely contained in the factor if, in which the particle dynamics 
enters only through the commutators. In the classical limit these commutators vanish and K becomes a simple 
numerical factor. Indeed in this classical limit one generally considers measurements of perfect precision, for which 
\f{^j)\^ ~^ H^j) ^'^d the factor K is unity. The expression (|4.2p then becomes the familiar form for classical Brownian 
motion (isj . On the other hand, in the quantum case, where the commutators do not vanish, it is clear that the 
particle dynamics is inextricably linked with measurement and, as a consequence, the symmetry property of classical 
stochastic processes does not hold. A second reason for the importance of this result is that it is convenient for 
calculation, as we shall illustrate in the example below. Before that, however, we give a brief derivation. 

Consider first the case n = 1. Using the expression (|3.29p in the definition (|4.ip of the characteristic function, we 
can write 

e(l) = ^y°° dx,r{xi - xih)}f{xi - x(ti)}e*'=i^^^ . (4.4) 

Making the change of variable xi xi + x{ti), we see that 

e(l) =if(l)(e*'=i"(*i)^, (4.5) 



where 



/OO 
dxir{xi)f{x,)e'''^-K (4.6) 
-00 

Recall that the sums in the expression (|4.3p are to be taken to be zero when j = n. 
Next consider 

^(1,2) = (r dx, r dx2/(l)V(2)V(2)/(l)e'('=^^^+'=^"^A 

- K{2) dx^rixi - x(ii)}e*{'=i-i+'=^-(*^)>/{xi - x{h)}J . (4.7) 

Here K{2) is exactly of the form (|4.6p but with the label "2" in place of "1 " . Remember also that with our shorthand 
notation the label also represents the instrumental parameters, so the measurement function / changes with the label. 
Next, we apply the generalized Baker-Campbell-Hausdorf formula (|B4p to write 



ei'^^^''^+^^'-'i^^)}f{xi - x{ti)} = f{xi - x{ti) + ik2[x{h), x{t2)]}e'i''''''+''''''^''^^ . 



(4.8) 
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With this and making the change of variable xi — > xi + x(ti) + k2[x(ti), x(t2)]/2i, we obtain 

^(1, 2) = K{1, 2) (^e^f^i^i^i) e^''2x{t2) ^kik2[x{ti),x(t2)]'^ ^ {'i-Q) 

where 

A (1,2)= / dxif {xi-k2 — }j{xi + k2 — }e 'A (2), (4.10) 

J-00 2i 2i 

which is of the form (|4.3p with rt = 2. As a final step, we use the Baker-Campbell-Hausdorf formula (|B3p to write 
gifeia;(ti)g»fc2:i;(t2)gfcife2[2;(ti),x(t2)] ^ gi{/cia:(ti)+fe2:E(t2)} ^nd obtaiu the form g3) with n = 2. 

For the general case, the argument goes in the same way. Assuming the form ( 14. 3p for smaller n, we write the 
definition (|4.ip in the form 

^(1, • • • , n) - K{2, . . . , n) - x{t,)}e'i'''-'+^'-- - x{h)}^ . (4.11) 

Then we bring the exponential factor to the right, using the theorem (14. 9p as in Eq. (|4.8p . Then, shifting the variable 
of integration and the using the Baker-Campbell-Hausdorf formula in the exponential factor, we get the form (j4.3p . 



A. Example: Wave packet spreading 

Here we consider the case of two successive measurements, each with a measurement function of the form (|3.16p . 
corresponding to a Gaussian slit. First, we consider a single measurement with 

With this we use the standard Gaussian integral (jBip . to evaluate the integral expression (|4.6p . We find 

K{l)^exp{^^alkl}. (4.13) 
Using the Gaussian property (jB6p to evaluate the expectation in (|4.5p . we find 

<Je^'=i"(*i)^ = exp{-i kl}. (4.14) 

With this, wc find 

C(l) = cxp{-ia2fc2}, (4.15) 

where we have introduced 

a'^{x^)+al (4.16) 
Finally, we invert the definition (I4.ip of the characteristic function to write 

/OC 11 

-lC(l)e-*^^ (4.17) 
-00 

This again is a standard Gaussian integral and we obtain the result 

W-(l) = -=L=exp{-|^}. (4.18) 

Thus the probability distribution associated with a single measurement is a Gaussian whose variance is the sum of 
that of the instrument and that of the underlying quantum state of the particle. We have presented the steps leading 
to this result in detail since these are the steps that will be used repeatedly in this and our later examples. 
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Consider now a pair of successive measurements, the first with measurement function of the form (|4.12p the second 
of the same form but with the index "1" replaced by "2" . Using the standard Gaussian integral to evaluate the 
integral in the expression (|4.10p for K{1, 2), we find 
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2) = exp{-ia?fc? - liai - t^^)kl}. (4.19) 

Then using the Gaussian property, we see that 

(^e'^''i^('i)+'''^it2)}'^ = exp{-i (x^) (kl + kl) - c{t2 - h)kik2}. (4.20) 

Here we have introduced the correlation (|A9p . Putting these together, using the expression (|4.2p for the characteristic 
function, we can write 

C(l, 2) = exp{-i((T2fc2 + 2aTpkik2 + r^k^), (4.21) 
where (note the misprint in the Eq. (7.18) of [l|) 



T 



1 ! 



[x{h),x{t2)f 



2 4a2 



arp = c{t2~ti). (4.22) 

Note that is the same quantity (|4.16p that appears in the single measurement function. The two measurement 
distribution function is given by 

W{1,2)= r ^ r ^^(i,2)e"'('=i"i+'=^"^). (4.23) 
J^oo 27r J^^ In 

With the form (|4.2ip of ^(1, 2) we can perform the integration using the multidimensional form (|B2p of the standard 
Gaussian integral. The result is 

^(M) . ^ ^„ exp{- -^-^-^;y;--%-;-^-^ }. (4.24) 

Here we remark first of all that the lack of symmetry of the quantum distribution is obvious: W(\, 2) ^ M^(2, 1). 
The exception is when the commutator vanishes. Note that the symmetry, or lack of it, is with respect to interchange 
of the labels, that is, one must interchange not only x\ ^ X2 and t\ ^ t2 but also a\ '^02- Another aspect of this 
asymmetry is that it is possible to make the last measurement one of perfect precision, that is, put 02 — 0. 

A second remark is that the time dependence is only through the time difference t2 — ^i. This is a general feature, 
independent of the form of the measurement function. It arises from the time-translation invariance of the equilibrium 
state. 

Finally, we remark that for widely separated times (t2 — i\ ^ 00) the correlation and the commutator for the 
oscillator vanish. Then we see that VF(1,2) W{\)W{2). This is a special case of the cluster property of the 
quantum joint distribution functions, a property they share with the classical functions. Whenever the time between 
any two successive measurements is large, the quantum joint distribution function factors into a product of distribution 
functions, the one corresponding to the earlier times, the other to the later times. The exception would be the case 
of a free particle, since in that case there is no approach to an equilibrium value of (a;^). 

V. THE PROBABILITY DISTRIBUTION 

In elementary quantum mechanics one interprets the absolute square of the wave function as the probability distri- 
bution of the particle position. That is, 'P{x, t) = \'ip{x, t)|^, where Vdx is the probability of finding the particle in the 
interval dx at time t. In our discussion this probability distribution becomes the conditional probability, specialized 
such that the second measurement is a perfect measurement. That, is, the probability distribution is given by 



Wll) 



0-2=0 
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Here we exhibit only the dependence on the final particle position and the time difference. (Due to the time-translation 
invariance of the equilibrium state only the time difference appears.) The picture we have is that the initial state 
is prepared by the first measurement, a measurement made on the equilibrium state, and that the second (perfect) 
measurement samples the state. Of course, in the special case of a particle not interacting with the bath and at 
temperature zero, this reduces to the elementary quantum mechanics prescription. 

If we use the expression p.29p for the joint probability distributions, we see that we can write 



While we can use the method described in Section [IV] to calculate W{\) and W{1, 2) and then put the results in the 
definition (|5.ip . it is just as well to evaluate the expression (|5.2p directly. For this purpose we introduce the integral 
expression for the delta-function. 



5{^) = / ^^"^'''^ (5-3) 



with which we can write, 



Now we use the formula (IB4p to write 

/(l)te— (*2)P/^j(l) = _ x,)f{x{h) - XI - GF)e-"(*^)^/^ (5.5) 

where we have used the relation (|A7p to write the commutator in terms of the green function, G = G{t2 — ti). Next 
we again use the integral expression (|5.3p for the delta-function to write 

r{x{h)-x,)f{x{h)-x,-GP) = / dx'J*ix[ + —)f{x[ 



2 ' 2 



2ttH 



Introducing this in (|5.5p and using the Baker-Campbell-Hausdorf formula (|B3p we can write 
(/(l)te-"(*^)^/V(l)) - dx[nx[ + ^)f{x[ ^) 



— OO 



'^^i(xt+x'^)P' /h (^^-i{x{tt)P' +x{t2)P) /h'^ ^ ^5 J-j 



The Gaussian property (jB6P allows us to write 

^^-.(.(U)P'4-.(*.))P)/.^ ^ exp{-Mi^^^±p±^}, 



(5.8) 



where c = c{t2 — ti) is the correlation (|A9p . Finally, with this result, the integral over P' in (|5.7|) is a standard 
Gaussian integral (|Bip and we obtain 

{{x^f-^)P^^ r , GP,^,^, GP^ 

(5.9) 



;(l)t,-..fe)PA;(l)^ = exp(-il^!L_^} |^ + ^)/(,; _ ^) 

exp{-(2l+£i)!-z^(a;i-fx;)£} 



Here we should recall that c — c{t2 — ii) is the correlation (|A9p and G — G{t2 — ii) is the Green function (|A5p . This 
is the key result of this section, valid for any form of the measurement function /. Dividing this result by its value for 
P = we get the integrand in the expression (|5.4p for the probability distribution. We next consider some examples 
corresponding to different choices of the measurement function. 
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A. Example: Wave packet spreading 

Again, we consider wave packet spreading with an initial measurement corresponding to a single Gaussian slit, the 
measurement function being that given in Eq. (I4.12p There is no need to evaluate the general expression (|5.9p since 
we already have the expressions (|4.18p for W{1) and (I4.24p for W{1, 2). Putting these in (|5.ip we can write 



This is a Gaussian distribution with center x{t) and variance uP'{t), where 



V{x,t) ^ }■ (5.10) 



X(t) — Xi 

a 

c{t) 



\xissUit)f c\t) 

-^""l 4^1 ^^-^^^ 

Here we we should again recall that c{t) is the correlation (IA9p . 

As a first consideration, we note that V(x, 0) is the probability distribution for the particle distribution immediately 
after the first measurement. Since c(0) = (^^) ^-nd the commutator vanishes at < = 0. we see that the center and 
variance of the initial distribution are 

m 



-^(0) ^ |y^. (5.12) 

As one can easily verify, this initial distribution corresponds to the product of a Gaussian distribution of variance 
(a;^) centered at the origin with one of variance centered at xi. That is, the initial distribution corresponds to the 
wave packet formed when the equilibrium state of the oscillator is passed through a Gaussian slit of width cti centered 
at xi. 



1. Free particle 

The free particle coupled to the bath in the absence of the oscillator potential corresponds to the limit (a;^) oo. 
The point here is simple: the oscillator force can be neglected near the center and the motion will be that of a free 
particle. Noting that c{t) — l^x"^) — s{t)/2, where s{t) is the mean square displacement and remains finite in the limit, 
we find that the center and variance (|5.1ip of the probability distribution become 

x{t) = X\ 

.^(t) ^.?+.(0-[^M!. (5.13) 

With the commutator expressed in terms of the Green function, this expression for free particle wave packet spreading 
corresponds to that obtained using path integral methods by Hakim and Ambegoakar [l6|. For a free particle not 
interacting with the bath and at temperature zero, in which case s[t) = and [a;(0), a;(<)] = iht/m, this reduces to 
the well known expression for wave packet spreading found in elementary quantum textbooks. 



2. Displaced ground state distribution 

Another limit of interest is that in which a\ — s- oo, while at the same time — > oo such that Xq = {x^) Xi/a\ is 
fixed. In this limit, the center and variance (15. lip become 

x[t) = y^a;o, 



'it) = {x'). (5.14) 
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The probability distribution (|5.12p therefore is that of a displaced equilibrium state. That is, 

r{x,t)^V,q{x~x{t)), (5.15) 

where 

is the probability distribution of the equilibrium state. Thus, the wave packet moves without spreading, the variance 
being that of the equilibrium state of the oscillator. The center is initially at xq and asymptotically approaches the 
origin as the equilibrium state is reached. We shall come back to discuss this result further in Section |V] when we 
discuss spreading of an initial coherent state. 

B. Example: "Schrodinger cat" state 

Here we consider the case where the initial measurement forms two separated wave packets. The first measurement 
function then has the form 

expl-fel^l + expl-feliyi!} 
f^^> [87ra2(l + e-'^V8-?)2]i/4 ' ^^■^^> 

With this form of the measurement function the integration in (|5.9p involves only the standard Gaussian integral 
(|Bip . Note that xi is the position of the center of the instrument, which should be chosen to be zero if we wish the 
wave packet pair to be symmetrically placed about the origin. We then find 



'(/(l)teH— (*^))^/V(l))\ _ exp{-^+*X2f} 



, Pd , {x^)(P . , GPd 
^ (cos — — + exp{ - „ p., „, ^ } cosh ^ „ 

where w'^ = 'uP{t2 — ti) is given in (|5.1ip and d = d{t2 — ii) with 



m . (5.19) 

While there is no difficulty evaluating the integral expression (|5.4p with this expression for the integrand, our interest 
will be in the limits of a free particle or, for the oscillator, a displaced ground state pair, in which case it is simpler 
to first evaluate the limits of the above expression and then evaluate the integral. 

1. Free particle 

As in the above example of wave packet spreading, we obtain the case of a free particle coupled to the bath in 
the absence of the oscillator potential by forming the limit (a;^) — > oo. Forming this limit of the expression (|5.18p 
then putting the result in the expression (j5.4p for the probability distribution and performing the integral with the 
standard Gaussian formula (jBljl we find 

1 fexp{-^^} , exp{-(^} 



V{x,t) = 



2 



+2/''-'-^'cosW°'-ff''|. (5.20) 

where now w'^{t) is given by the free particle form (j5.13p and we have used the relation (jA7p to reintroduce the 
commutator. In this expression a{t) is the attenuation coefficient, given by 

a(0=exp{-^-Ji^}. (5.21) 
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This expression for the probabiUty distribution is the same as that derived in an earher brief communication fl7| . 

This probabihty distribution is the sum of three contributions, corresponding to the three terms within the braces. 
The first two are probabihty distributions of the form (|5.12p corresponding to a pair of single shts positioned at ±d/2, 
while the third term (that involving the cosine) is an interference term. The attenuation coefficient is a measure of 
the size of the interference term and is defined as the ratio of the amplitude of the interference term to twice the 
geometric mean of the other two terms. The point here is perhaps best seen if we look first at the case of a particle 
without dissipation and at zero temperature. Then s{t) = and [a:(0), a;(i)] = iht/m. The resulting probability 
distributions are shown in Figure ??. There V{x,0) is the initial probability distribution, which depends only on the 
initial measurement function and is therefore the same whether or not there is dissipation. In this initial distribution 
the interference term corresponds to a miniscule peak at the origin, so small that it does not show in the plot. In the 
same figure V{x,t) is the distribution at a time t such that the width of the individual wave packets has increased 
by a factor of roughly 3, while Vo{x,t) is the same distribution but with the attenuation factor a(t) set equal to 
zero. We emphasize that at this later time the amplitude of the interference term is of the order of that of the other 
terms, despite the fact that initially it is negligibly small. The difference between 7'(x, i), where the attenuation factor 
is unity and the interference term is present, and Vo{x,t), where the interference term is absent, is what is called 
decoherence. Thus, the attenuation coefficient corresponds to the traditional measure of coherence in terms of the 
relative amplitude of an observed interference term.[l8| 

In the presence of dissipation the attenuation factor decays rapidly when the separation d of the wave packets is 
large. Here by "large" we mean not only large compared with the slit width cti but also large compared with the mean 
de Broglie wavelength A — hjra^J {x^). To obtain the short time behavior in this case, we note that, as we have seen 
in Section. IlIIl for short times the mean square displacement s{t) = (^i^) t^. Then, since for short times w'^{t) = af, 
we see that for short times, 

a(t) ^ e-*'/2"S (5.22) 

where Td, the decoherence. time, is given by 

In the high temperature case, where (i^) = kT/m, this is the result for the decoherence. time obtained previously 
[l3|, but the result holds equally well at zero temperature, where (i^) is given in Eq. (|A29[) . In either case, the 
decoherence. time is very short when the separation d of the pair of wave packets is large. 



2. Displaced ground state pair 



As in the above example of wave packet spreading, we obtain a relatively simple expression for the oscillator case 
in the limit cr^ — > oo and d oo such that do = (x^) d/af) is fixed. Forming this limit of the result (|5.18p and then 
evaluating the expression ( 15. 4[) for the probability distribution, we find 

- 2(l + e-W)) {^-'^(^-^)+^--^- + ^) 

2ait)e-'-i*y'i^')v^,ix)cos , (5.24) 

4i (a;^) J 



where Vcq is the equilibrium distribution (|5.16|) and now 



d{t) = ^do, (5.25) 
\x ) 



while here the attenuation coefficient is given by 



In connection with this result we remark first that initially the distribution is of the same form as in the free particle 
case, with two Gaussian peaks at ±c?o/2 and a miniscule central peak. The difference is in the motion: the two peaks 
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drift back and forth against each other without spreading, eventually arriving at the origin. The interference term 
therefore has a different effect (sometimes called a "quantum carpet" [l^) but we have nevertheless introduced the 
attenuation coefhcient in the same way as the ratio of the coefficient of the cosine to twice the geometric mean of the 
first two terms. For very short times, c{t) = (x^) — ^ ^i^) and [x{0),x{t)] ~ iht/m, so a(t) is of the same form 
(|5.22[) but now with 

Td = ^ ^ 5.27 

Note that the uncertainty principle tells us that (i^) (x^) > ?i^/4, so the argument of the square root is always 
positive.. 

VI. THE WIGNER FUNCTION 

The Wigner function is the analog of the probability distribution in which the second measurement is a perfect 
measurement of position and momentum. Of course, a quantum particle cannot have simultaneously a precise position 
and momentum, so this last cannot be a proper quantum measurement of the general form (j3.7p . The Wigner function 
is therefore not a probability distribution but rather what is called a "quasiprobability distribution function" . To get 
the Wigner function corresponding to the density matrix at time t2, which we denote by W{q,p', t2 — ii), we make the 
replacement /^(2)/(2) "(^(g — x{t2)S{p — 'mx{t2))" — F{q — x{t2),p— mx{t2)), where 

F(,,,).^y_dgy_dPe^(--^^)/^ (6.1) 



That is, in place of the general formula (15. 2p for the probability distribution wc have its generalization to quantum 
phase space given by the general formula 

, , , {fHl)F{q~x{t2),p-mx{t2))fil)) 

m,.P-M~t.) - ^T(T)7(T)^ 

dQ p dP^,(p,+Qrt/;,(/ni)e-^(^(*^)^+'"*(*^)'3)A/(l)) 



2nhJ_^2nh {fHl)f{l)} ' ^'^ 

Note first of all that as a classical function F{q,p) = 6{q)S{p) . This expression would be unsatisfactory for our 
purposes, since in the general formula (j6.2p the arguments of the delta functions would not commute. The integral 
expression ()6.ip corresponds to the Fourier-von Neumann representation of the classical operator. Second, we note 
that the momentum operator is interpreted as the mechanical momentum mx . This, as we have noted above, is in 
accord with the macroscopic description of a dissipative system. The point here is that the canonical momentum is an 
operator of the microscopic description, which for the same macroscopic description may or may not be equal to the 
mechanical momentum [20|. Finally, we should emphasize that the Wigner function is not a probability distribution 
as discussed in Sec. IIIII For this reason we use a calligraphic W to help keep this in mind. This formula for the 
Wigner function is unique in the sense that it satisfies certain general requirements such as that it be a real function, 
that the integral over p or q must give the corresponding probability distribution in position or momentum, etc. For 
a thorough discussion of the Wigner function as it has appeared in the literature, see the review article of Hillery et 
al. [m (see especially their Eq. (2.45). 

It is convenient to introduce the Fourier transform of the Wigner function, 

/>00 /"CxD 

yV{Q,P;t)^ dq dpe-'^P''+'^P^'''W{q,p-t). (6.3) 



This Fourier transform is what in the literature is called the "characteristic function" [21\ . We adopt this convenient 
terminology, but warn that this Wigner characteristic function should not be confused with the quantum analog of 
the characteristic functions of classical probability introduced in Section [IVI The inverse Fourier transform is 



W(g,p;i) ^ f ^^y^i^Q^p.t)e^(P^+Qv)/\ (6.4) 

Comparing this with the general formula (|6.2p . we obtain by inspection a simple formula for the Wigner characteristic 
function: 

W{Q,P-M^h) = (;t(i)/(i)) ■ (6.5) 
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This is the key result of this section, vaUd for any form of the measurement function /. As we next shall show, it 
allows us to readily calculate the Wigner function for a variety of examples. 

A. Example: Equilibrium Wigner function 

As a first simple example we consider the equilibrium Wigner function, which we denote by yVoq('Z,p), and which 
we get when we make no initial measurement. This corresponds to /(I) — > 1, and for this case the formula (|6.5p for 
the Wigner characteristic function becomes 

VVoq(Q, P) - ^e-*(^(*^)^+™*(*^)'3)/'"'^ . (6.6) 

Using the Gaussian formula (jB6|) we find 

Weq(g, P) = exp{-^((a;2) + (x^) Q')}, (6.7) 

where we have used the fact that {xx + xx) — 0. With this in (16. 4p . the integrals are standard Gaussian integrals and 
we find for the equilibrium Wigner function, 

Here {x^) and (i^) are given in the expressions (jA30|) and (|A29I ). The familiar weak coupling form, well known as 
the equilibrium solution of the master equation, results if we recall the relations for weak coupling given in Eq. (IA39|1 . 

B. Example: Motion of a coherent state 

Coherent states are generally defined for the free oscillator by operating on the oscillator ground state with the 
general displacement operator eyi'p{i{mv{)X — xop)/h}. The resulting coherent state corresponds to a displaced ground 
state, centered at xq and moving with velocity vq . Here we define a generalized coherent state for an oscillator 
interacting with a linear passive heat bath. The corresponding density matrix is obtained by acting on the equilibrium 
density matrix with the measurement function 

/(I) = exp{im(uox(ii) - xox{ti))/H}. (6.9) 

With this the expression (|6.5p for the Wigner characteristic function becomes 

yV{Q, P; t2 — ti) ~ ^^-im{vox(ti)-Xox(ti)/h^-i{x{t2}P+mx{t2}Q)/h^im{vox{ti)-Xox(ti)}/h^ ^ (6.10) 

Use the Baker-Campbell-Hausdorf formula (|B3[) to write 

^-im(vox{ti)-xox{ti))/h^~-i{x{t2)P+mx{t2)Q)/h^im{vox{ti)-xox{ti))/h 
— ^-i{x{t2)P+mx(t2)Q) / h ^mlx{t2)P+mx(t2)Q .vox{ti)-xox{ti)]/ j-g ^2) 

With this, we find 



W(Q, P; t) = exp{-^i (.x(t)P + mt;(i)Q)}Weq(Q, P), (6.12) 



where, using Eq. (|A7p to express the commutator in terms of the Green function, 

x{t) — ■mG{t)xQ + ■mG{t)vQ, 
dx(t) 

v{t) = ^ 7nG{t)xa + mG{t)va. (6.13) 

at 

The Wigner function is given by the inverse Fourier transform (16. 4p . We find 

W(g,p;i) = Wcq(q-x(t),p-TOv(t)). (6.14) 
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Thus the Wigner function corresponding to an initial coherent state has the form of an equihbrium Wigner function 
whose center moves according to the equations (|6.13p . 

We should point out that the motion (|6.13p of this center is not the solution of the mean of the quantum Langevin 
equation (jAip . Rather, it is the solution of the mean of the initial value Langevin equation [23 |. 

+ f dt'Kt - t')^^ +Kx = -fi{t)x{0), (6.15) 
di'^ Jo at' 

with initial data x{0) = xq and v{t) = vq. The effect of the term on the right hand side can be seen clearly in 
the Ohmic case, where the Green function has the form (jA20|) . With that form we see that x{t + 0~^) = xq, while 
v{t + 0~^) — Vq 7x0. Thus the center of initial distribution in the qp plane makes a jump in the p direction, down 
(up) if xo is positive (negative), after which the motion of the center is that of a damped harmonic oscillator. At all 
times the shape of the distribution is that of a displaced thermal equilibrium state of the oscillator. 
The probability distribution is obtained by integrating over p, 

/OG 
dpW{x,p;t). (6.16) 
-00 

With the expression (|6.14p for the Wigner function in which Wcq is of the form (|6.8p this becomes 

r{x,;t)=req{x-x{t)), (6.17) 

where Vcq is the equilibrium distribution (|5.16l ) . This is exactly the form we encountered above in the example of 
wave packet spreading, where the parameters for the displaced ground state are given in Eq. (|5.14p . In either case 
the probability distribution is that of a displaced ground state. The difference lies in the motion of the center, which 
is temperature independent in the case of the coherent state but has a temperature dependent form for the displaced 
ground state. In either case, of course, the distribution approaches that of equilibrium for long times. The lesson we 
learn here is that the time dependence of the approach to equilibrium depends on how the initial state is formed. In 
Figure ?? we plot x{t) for the two cases, the displaced ground state motion being calculated at zero temperature. 
The parameters chosen were "f/tOQ — 10/13 and Q/ujq — 5, but despite this rather strong coupling, there is not much 
difference between the two curves. 



C. Example: Coherent state pair 

The idea here is to form an initial state like the "Schr odinger cat" state discussed in Sec lVBl There the initial 
state was prepared with a pair of Gaussian slits. Here we consider instead an initial state which is a superposition of 
two separated coherent states. This is accomplished with a measurement function of the simple form 

/(l) = cos!^, (6.18) 

which results in a superposed pair of generalized coherent states, centered at a; = ±c?/2 and each with zero velocity. 
The expression (j6.5p for the Fourier transform of the normalized Wigner function at time t2 therefore becomes 

/cos in:^£ihlp-i(x(t2)P+mx{t2)Q)/fi „„g mdxjti) \ 

W{Q,P-M-ti) = ^ — — ^ 6.19 

Now, using again the Baker-Campbell-Hausdorf formula (|B3p . we see that 

mdi(ti) ^_,(,(,^)P+„i(,^)Q)/;, mdxjti) 

- ^ (o^-i(x(t2)P+mx(t2)Q)/h md{GP + mGQ) 

_^^-i(x(t2)P+mx(t2)Q—mdx(ti))/h _|_ ^—i{x{t2)P+mx(t2)Q+mdx{ti))/h^ j-g 
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where, to shorten the expression, we have used the expression (|A7p to express the commutator in terms of the Green 
function. Putting this in (|6.19p and using the Gaussian property (|B6[) we find 

~ cxp{-((x^)p2 + m^(i:^)Q^)/2nn / md{GP + m,GQ) 
= l + exp{-mMi:^)rfV2^^ 1,'°' 

^ m?(P{x^)^ md(sP + msQ)\ 
+ cxp{ ^} cosh ^^ J , (6.21) 

where G = G(t) is the Green function (|A5p and s = s{t) is the mean square displacement (|A8p . Putting this in 
the integral is a two dimensional standard Gaussian (jB2p and we find 



1 I , mdd m'^Gd, 

= 2(1 + exp{-m^ (A^) d72^n) T"^^'^ ' ' 

+Wcq(g + ^,P+^— ) 

+2e-^(*) Weq(q, p) cos $(q, p; t) } , (6.22) 
where Weq('?,p) is the equilibrium Wigner function (j6.8p and we have introduced 



A(t) = 



2^2 I 4(.T2)(i2) 4 (±2)2 



*<-')^(i5f+(Sy)4 (-3) 

This expression for the Wigner function is identical with that obtained by Romero and Paz using path integral 
methods [23]. 

Viewed in the qp plane, the expression (|6.22p for the Wigner function shows three peaks: an outlying pair in the 
form of single coherent states, centered initially at g = |, p = 0, and an interference peak, centered at the origin and 
modulated by the factor cos$. Initially, since A{Q) = 0, the amplitude of the interference peak is twice that of either 
of the two outlying peaks. In general the interest is in the case of a widely separated coherent state pair, that is, 
d ^ (a;2). In that case for very short times A{t) becomes large and the interference peak is practically zero. This 
disappearance of the interference peak is the phenomenon of decoherence. as seen with the Wigner function. To be 
more explicit, with the expansion (|A26p in the expression (|A24p for s{t), we can readily evaluate this expression for 
short times. In the limit of small bath relaxation time this takes the simple form: 



A{t) = -^^ogl, (6.24) 




where 

rd = rW^. (6.25) 

Here C is the friction constant and r is the bath relaxation time, the parameters in the single relaxation time model. We 
can interpret Td, the time of the order of which the central peak in the Wigner distribution vanishes, as a decoherence. 
time. Note that this expression for the decoherence. time is very different from that for the displaced ground state 
pair given in Eq. (|5.27p . despite the similarity of the initial states. The qualitative nature of the phenomenon is the 
same: the rapid disappearance of an interference term. 

Since the Wigner function is not directly observable, we should consider the probability distribution. That is, we 
put the expression (j6.22p for the Wigner function in the integral (|6.16p for the probability distribution to obtain 



2(l + exp{-J(.x2)d2/2;,2|) (^oq(. - mG^) + V.,{. + mG^) 



I rn^G'^d'^ 1 ^ / \ msdx \ ,^ 
2aexp V.,[.) cos , (6.26) 
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where a = a(t) is the attenuation coefRcient, now given by 

As in our discussion of the Schrodinger cat state in Sec. |V] , the attenuation coefficient is defined as the ratio of the 
coefficient of the cosine term to the geometric mean of the first two terms and corresponds to the traditional measure 
of coherence. Its disappearance is the phenomenon of decoherence. as seen in the probabihty distribution. But here 
the initial value of the attenuation coefficient, 

a(0) = exp|-(^^^-lj^|, (6.28) 

is already vanishingiy small for large separations (the uncertainty principle tells us that the factor in the exponent 
is necessarily positive). Indeed, for the Ohmic model, the mean square velocity is logarithmically divergent, so the 
attenuation coefficient is identically zero for all times. In any event, we would say that as seen in the probability 
distribution, the decoherence. for a coherent state pair occurs initially and there is no notion of a decoherence. time. 

The earliest discussion of this problem of a displaced pair of coherent states at zero temperature was that of Walls 
and Milburn,[23| who based their discussion on the master equation. As we note in the last paragraph of appendix 
El this would correspond to the Weisskpof-Wigner approximation. In this approximation 

^W = ^(l-e-*) (6.29) 

and 

a{t) = exp{-yl(i)}. (6.30) 

Note first of all that in this Weisskopf-Wigner approximation the decay of coherence is identical whether viewed in 
the Wigner function or in the probability distribution. However, the short time behavior is very different from that 
of the exact result. As an illustration, in Figure ?? we compare the short time behavior of the exact result with 
that of the Weisskopf-Wigner approximation as well as that of the weak coupling approximation. In this range the 
weak coupling expression for A{t) is just twice that from the Weisskopf-Wigner approximation, while both are much 
larger than the exact result and would give a correspondingly much shorter estimate of the decoherence. time. The 
conclusion to be drawn is that the master equation can give misleading results at short times. 



D. Example: Squeezed state 

The squeezed state that appears in the quantum optics literature is obtained by operating on the ground state of 
the free oscillator with the squeeze operator [25], 

S = exp{^{e~'''a^-e''a^^)}, (6.31) 



where r and 9 are real parameters and a ~ {uiujqx + ip) / ^/2mhujo is the usual annihilation operator for the free 
oscillator. Here we consider the so called ideal squeeze operator, corresponding to = 0, and again replace the 
canonical momentum with the mechanical momentum. The ideal squeeze operation would therefore correspond to an 
initial measurement operator of the form 

777 T 

/(I) = exp{i — {x{h)x{h) + x{h)x{h))}. (6.32) 
Since this is a unitary operator, we see that 

^/t(l)e"*^''^*"^^+"*^*"^'^^/^/(l)^ = (^e"*('^('^'*")^+™'*('''*")'3)/'*^ , (6.33) 
where we have introduced the operator 

X(r;i2)=/^(l)x(t2)/(l). (6.34) 
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To evaluate this operator, form the derivative with respect to r, 

^ -mG{t2~h)X{r;ti)+mG{t2-h)X{r;ti), (6.35) 
where we have used the relation (IA7p to express the commutator in terms of the Green function.. Also, we have 

^^^pM. = -m,G{h - h)Xir; t,) + mG{h - h)X{r- h). (6.36) 
or 

If we set ^2 = ti and use the fact that G(0) = and G(0) = l/m, we find from (|6.35p that 

X(r;ii) -a:(^l)e-^ (6.37) 

while from (|6.36p we find that 

X{r;ti) ^ i(ti)e\ (6.38) 

Here we should emphasize that we use the single relaxation time model, for which G(0) = 0. For the Ohmic model 
there would be an extra term. Putting these results in (|6.35p and (j6.36p and integrating, we find 

X{r;t2) = x{t2) - {I - e''')mGx{ti) + [e'' - l)mGx{ti) 

X{r;t2) = x{t2)-{l-e-'')mGx{h) + {e'' -l)mGx{h). (6.39) 

where G = G(t2-ii). 

Since X{r; 12) and X{r] ^2) are linear in the operators x(ti) and x{t-C), they have the Gaussian property and we can 
use the identity (|B6|) to evaluate the expression l|6.33p . With this result the Wigner characteristic function (|6.5p can 
be written in the form 

Win Pi \ s AiiP^ + 2Ai2QP + A22Q\ ... 
M/(Q,P;t2) = exp{ — }, (6.40) 

where for this present example, 

An = {X^r-M)) 

= [1 - (1 - e--)mGf {x') + {e- ~ ifm^G^ (x^) 

+ (1 - e-'')TOGs + (e"" - 1)toGs, (6.41a) 
A12 = - (^X{r;t2)X{r;t2) + X{r;t2)X{r;t2)) 
= {l~e-'')mG]{l-e-~'')m^G{x^) 

+ {e'- - lym^GG (x^) + ^(1 - e-'')™^ (gs + Gs) 

+ ^{e'' - l)m^ (Gs + Gs) , (6.41b) 

A22 = m^(^X^r:t2)] 

= [1 + (e-^ - l)^m^G^]m^ (i^) + (1 - e-^)^m^G^ (x^) 

. + (€'■- l)m^Gs + (1 - e-'')m^Gs. (6.41c) 

Here we should again recall that s — s(<2 — ^i) is the mean square displacement (jA8| and G — G{t2 —ti) is the Green 
function (KE\i . 

Forming the corresponding Wigner function, we find for the squeezed state, 

,AU ^^ 1 r AnP^ - 2Ai2Pq + A22q^ 

271^/^11^22 - A?2 2(^11^22 - Ai2) 

In Figure. ?? we plot constant density contour for this function in the plane of the dimensionless variables u — q/ (a;^) 
and V — p/m(^x^y The dashed circle corresponds to the equilibrium state, the state just before the initial squeeze 
as well as the state at long times. The contour marked (0) corresponds to the initial squeezed state. In the course 
of time this contour rotates clockwise. The contour marked (1/4) is that corresponding to a quarter period, while 
that marked (1/2) is that corresponding to a half period when the squeezing is much reduced. The relatively strong 
coupling chosen, j/ujo = 10/13, emphasizes the effect: dissipation leads to a loss of squeezing. 
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E. Example: "Schrodinger cat" state 

Here we consider the Wigner function when the initial measurement corresponds to the measurement function 
(j5.17p . The calculation goes exactly the same as that beginning with (|5.5p in the previous section, so we shall simply 
quote the results. For general / we find 

1 , cP + mcQ 

Here c = c{t2 — ti) is the correlation (jA9[) and G = G(t2 — ii) is the Green function (jA5[) . With the form (|5.17p for 
the measurement function, this becomes (remember we must choose xi = if we wish the initial wave packet pair to 
be placed symmetrically about the origin). 



/ xi=0 

AiiP^+2Ai2PQ+A22 



(l + e-<iV8-?)^27r((a;2)+a2) V 2((a;2) + a^)?! 



+e-Vs.?eoshi^^±^l, 
where in this present example 



4al 



(6.44) 



An - (.^) + 



4^2 (2;2) 
'h^GG cc 



Acrf (a;2) + erf 



^22 = ( (i') + ^ - , / . 1 ■ (6.45) 



4(7^ 



The Wigner characteristic function is therefore 



l+exp-[ 8cr^((^2)+<^2)J- I. V\ / T 

+ ^-P^- 8.2(V)+.2) }cosh^ ^-^1 . (6.46) 

With this form of the Wigner characteristic function there is no difficulty evaluating the inverse transform (16. 4p to 
obtain the corresponding Wigner function. However, the interest will be in the limits of a free particle or, for the 
oscillator, a displaced ground state pair, in which case it is simpler to first evaluate the limits of the above expression 
and then evaluate the inverse transform. 



1. Free particle 



We obtain the case of a free particle coupled to the bath in the absence of an oscillator potential by forming the 
limit {x^} ^ oo. In forming this limit we should recall that c(t) = {x"^} ~s{t)/2, where the mean square displacement 
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s{t) is finite in the free particle limit. We find 

cxp{- 

l + exp{-^} 



„„„r A^^P'+2A^2PQ+A22Q' ^ ( p . 



, (P ^ , iGP + mGQ)d\ 
+ exp{-g^}cosli^ )■ (6-47) 



In this free particle case, the expressions (|6.45p become 



All = crf + s 



4a2 ' 



/ s ti^GG 

Al2 = "I I - + — 2- 



2 Aoi 



Forming the inverse Fourier transform (|6.4p we obtain 



A22 = m2((i') + ^V (6.48) 



2(i + e-'iV8-f) 



2'"' ' 2' 



+2 exp{-A(t)}Wo(q,p; t)cos$(q,p; t)} , (6.49) 
where Wo is the Wigner function corresponding an initial measurement forming a single wave packet at the origin, 

r _ A22g^-2Ai2<?p+Aiip^ 1 

Wo(g,p;t2)- — . ^ ^ • (6.50) 

27r^AiiA22 - Al^ 

In this example the phase $ is given by 

. . (GA22 - mGAi2)g + (mGAii - GAi2)p M 

"^'^'-^'^ = A11A.2-AL 4^ (^-^^^ 

and the quantity A by 

AW ^ ' 7^ ; ' i;. (6.52) 

^ ' A11A22 - ^ ' 

As in the case of the coherent state pair, The Wigner function for the free particle "Schrodinger cat" state shows 
three peaks, an outlying pair centered at q = ±d/2,p — and an interference peak centered at the origin. However, 
in this free particle case 

A{0) = -Ai^^JL ^ f > (6.53) 

where A is the mean de Broglie wavelength in equilibrium, 

A = (6.54) 

Therefore, the amplitude of the interference peak will initially be vanishingly small whenever the separation d of the 
wave packets is large compared with both the slit width ai and the mean de Broglie wavelength A = hjm^J (i^). 
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2. Displaced ground state pair 



Again we form the limit ct^ oo and d —> cxd such that do = {x^) d/o'i is fixed. The coefficients (|6.45p become 

An = (x2), Ai2 = 0, A22^m^x^).. (6.55) 
The Wigner characteristic function (|6.46p then becomes 



cxp{-— 2h-^ } f {cP + mcQ)da 

cos ■ 



WiQ,P;t) 



l + exp{-d2/8(x2)} \™ 2{x^)h 
[GP + mGQ)do 



With this, the Wigner function is 
W{q,p;t) 



3xp{ -dl/8(x'^)} cosh ■ 



4(2;2) 



2(l + cxp{-^} 



cdn 



2(x2 



mcdo 



+'^cq(<Z + ^^,P+ J^) + 2e ^>Veq(<z,p)cos$ 



2 a;2 



where yVcqiq,p) is the equihbrium Wigner function (|6.8p and 




Again, as in the free particle case the initial value of A is not zero. 



AiO) 



4 



(X2) 



4m2 (±2) (a;2) 



> 0, 



(6.56) 



(6.57) 



(6.58) 



(6.59) 



and the interference term is vanishingly small. 

This situation, in which the initial state is a "Schro dinger Cat" state formed by passing the particle through a 
pair of Gaussian slits, is to be contrasted with that described in Sec. IVI CI in which the initial state is prepared by 
displacing the equilibrium state to form a coherent state pair. The Wigner functions, given by Eqs. (|6.58p and (|6.22p . 
respectively, are identical in form, but in the "Schrodinger Cat" case the interference peak is initially vanishingly 
small and remains so for all time, while in that of the coherent state pair the interference peak is initially twice as 
high as the outlying peaks, becoming vanishingly small only after a short relaxation time. The reverse is true for 
the probability distributions, given in Eqs. (j5.24p and (j6.26p . respectively. That is, the interference term in the 
probability distribution is vanishingly small at all times for the coherent state pair, while for the "Schrodinger Cat" 
state the attenuation coefficient multiplying the interference term vanishes only after a short decoherence. time. We 
must conclude that the notion of decoherence. time is arbitrary, depending on the situation and how one chooses to 
view it. 



VII. CONCLUDING REMARKS 



The quantum probability distributions are measured distributions. That is, they depend explicitly on the parameters 
of the measurements. This is seen clearly in the general formula (j3.29p for the n-point distribution, where /(j) is the 
measurement function for the j'th measurement. This is also seen in the in the case of quantum Brownian motion in 
the expression (|4.2p for the characteristic function. There it is seen that the characteristic function can be factored, 
with a factor K{1, ■ ■ ■ , n) multiplying a quantum expectation independent of measurement. The factor K{1, ■ ■ ■ ,n) 
contains the measurement parameters and depends upon the dynamics through the non-equal-time commutator. 
In the classical limit, where this commutator vanishes, this factor becomes a numerical factor independent of the 
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dynamics that can in practice be taken to be unity. The result is the famihar expression for the classical characteristic 
function, [l^ 

An important part of this work has been to demonstrate that the general expression (|4.2|) for the characteristic 
function can be very useful for practical calculations. The use of this formula together with its specializations (j5.9p 
to the probability distribution and (j6.5D to the Wigner function has been illustrated with a number of examples, each 
of which is an important application. Among the results we point out the expression (I5.13|) for wave packet spreading 
in the presence of dissipation, a generalization of the well known expression found in elementary textbooks. Another 
result is illustrated in Fig. ?? where the disappearance of squeezing in the presence of dissipation is illustrated. 
Finally, a comparison of a "Schr odinger Cat" state formed by passing the particle through a pair of Gaussian slits 
with the nearly identical state formed with a coherent state pair shows that a quantitative measure of decoherence. 
depends on how the state is formed. 

In Sec. IVII we discuss the Wigner function. Some authors prefer instead the density matrix element in the coordinate 
representation, which is given by a kind of half-Fourier transform: 

/oo I I 

dpW{^^,p;t)e^^^-^'^P/\ (7.1) 

Clearly, the Wigner function and the density matrix element contain the same information. Indeed, it is not difficult 
to see that the central interference peak in, say, the Wigner function (|6.49p corresponding to a coherent state pair 
becomes a pair of off-diagonal peaks in the density matrix element. We prefer the language of the Wigner function 
since it is always a real function that in the classical limit becomes a real probability distribution. 

The examples are all in one dimension. The reader should be aware that this is not a necessary restriction, but has 
been made to keep the discussion within bounds. The generalization to higher dimensions is straightforward. All that 
one must keep in mind is that the fluctuating force operator is independent in the different directions and that the 
effects of dissipation (the fluctuating force and the memory force) are independent of the applied force, [l^l Finally, 
we have restricted the discussion to the single relaxation time model of dissipation and its limiting Ohmic case. Again 
this has been done to keep the discussion within bounds. There is no problem with the discussion for more general 
models such as the coupling to the blackbody radiation field. [20j 



APPENDIX A: QUANTUM BROWNIAN MOTION 

1. Quantum Langevin equation 

Quantum Brownian motion for an oscillator coupled to a heat bath at temperature T is described by the quantum 
Langevin equation, 

mx+ f dt'^i{t-t')x{t') + Kx ^ F{t). (Al) 

J — oo 

This is a Heisenberg equation for the position operator x{t). On the right hand side F{t) is a Gaussian random 
operator force, with mean zero, {F{t)) — 0, and with autocorrelation and commutator given by 

l(Fit)Fit')+F{t')Fit)) 

In these expressions fl is the Fourier transform of the memory function, 

/i(z) = / dt^i{t)e'^\ Im{z} > 0. (A3) 

"'0 

It is a consequence of the second law of thermodynamics that jl{z) must be what is called a positive real function: 
analytic with real part positive everywhere in the upper half plane. 

In our present discussion we take the view that the above is a macroscopic description, which is complete as it 
stands. For a thorough discussion, including the derivation from a number of microscopic models, we refer to a paper 
of Ford, Lewis and O'Connell [20j . 



^ (iwRe{/i(w -I- iQ^)}uj coth — — cosLjj{t — t'), 

2kT 



IT 

h 

in 



du}Re{p.{uj + iQ+)}LJsmLj{t-t'). (A2) 
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The solution of the quantum Langevin equation (|Aip can be written 

x{t) ^ I dt'G{t~t')F(t'), (A4) 

where the Green function G{t) is given by 

1 /"^ 

G{t)::^— duja{uj + iQ+)e-''^\ (A5) 

in which a{z), the response function, is given by 

oi{z) = ^ ^ . (A6) 

With this solution, we can obtain the following expressions for the commutator, 

2h 1"°° 

\x{t),x{t')] = — / dujlm{a{uj + iO+)}s\nuj{t-t') 

= ih{G{t' -t)-G{t-t')}. (A7) 

The Green function vanishes for negative times, while G(0) — and G{Q) = 1/m. We see therefore, that the canonical 
commutator, [x,p\ = ih, holds with p = mi, the mechanical momentum. (The canonical momentum may or may not 
be the same as the mechanical momentum, depending on the form of the microscopic Hamiltonian.) 
Also of interest is the mean square displacement, 

Sit-t') EE {ixit)-xit')r) 

= 2{x^)-2c{t-t'). (A8) 

Here c{t) is the correlation 

c{t) = ^{x{t)x{0)+x{0)x{t)). (A9) 

Using the solution (|A4p or, more directly without recourse to the Langevin equation, the fluctuation-dissipation 
theorem of Callen and Welton [2^ , we obtain the following expression for the correlation 

c{t) ^ - I dcjlm{a(cj + i0+)}coth— ^^coswt. (AlO) 
TT Jo 2fcT 

With this, we also have 

sit)^ — / dwlm{a(w + i0+)}coth (1-coswt). (All) 

TT 7o 2fcr 

Here we should note first of all that for very long times the correlation vanishes and 

s{t)^2{x^). (A12) 

The exception is the free particle case, where K — Q and consequently (x^) = oo. In that case s{t) grows for long 
times without limit, with a time dependence that depends on the model as well as the temperature [23, HI]. On the 
other hand, for very short times we can expand the cosine to obtain the general result 

s{t)~{x'')t'. (A13) 



2. Explicit expressions 



Here we obtain explicit, closed form expressions for the Green function G{t) and the mean square displacement 
s{t). For this purpose, the model of choice for most applications, due to its simplicity, is the Ohmic model, for which 
/i(z) = the friction constant. While it is adequate for most purposes, this Ohmic model is singular, particularly 
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at short times or high frequencies, so we here consider a more general model in which the higher frequencies are 
suppressed. The simplest of these is the single relaxation time model, for which 

f,{z) = (A14) 

1 — IZT 

Here t is the relaxation time (at times called the bath correlation time) which we assume is small in the sense that 
Cjjm <C 1. Putting this form of /x in the response function (jA6[) and replacing the parameters r, ^ and K with 
parameters , 7 and through the relations 

C = m7 \^ \ K^mujl——, (A15) 



S2 + 7 (S2 + 7)"^ i2 + 7 

we obtain the convenient form 

a[z) = — , ^ ^. (A16) 

Note that — ifi is the pole of the response function far down on the negative imaginary axis. For small relaxation time 
57 ~ 1/r we have the expansion 



T m 77?/ ~ 777'' 777-^ 

,2 



The Ohmic model, for which t 0, therefore corresponds to 17 ^ 00, while 7 Qjm and Wg K/m. 



a. Green function 

With the form (jA16P of the response function, we evaluate the integral in the expression (|A5P for the Green function 
by deforming the path of integration into the lower half plane, picking up the residues at the poles of the response 
function. The result is 



Vfi — •ffl + 777.ti;i 



g-7t/^:_^::^^ (AI8) 



where 



^^l - 7V4. (A19) 



In the limit, 17 — > 00, this becomes the familiar Ohmic Green function 



<-^Ohmic I — e . 



^Ohmic(i) = e'^'/^ (A20) 

777^1 

For reasonable choices of the parameters the Ohmic Green function is very little different from that for the single 
relaxation time model. The difference between the two models is more apparent in the second derivative of the Green 
function, which wc show in Figure ?? for the two models. The difference is small except for short times, where G(0) 
vanishes, while Gohmic(O) — —C^jw? is finite. 



6. Mean square displacement at high temperature 

First we consider the mean square displacement in the high temperature limit, in which in the expression (jAlOp we 
replace the hyperbolic cotangent by the reciprocal of its argument. If we compare the resulting expression with the 
expression (|A5|) for the Green function, we see that in this high temperature limit, 

s{t) = 2kT I dt'G{t'). (A21) 
Jo 

The long and short time behavior of the mean square displacement are of the general form given in Eqs. Cn2| and 
(jAlSp . Note that in this high temperature limit {x^) and (i^) are given by the classical equipartition formulas. 
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c. Mean square displacement at zero temperature 

We consider now the effect of zero-point oscillations on the mean square displacement. At temperature zero, we 
replace the hyperbolic cotangent in the expression (jAlOp by unity. Then if we write, 



we can write the expression (jAlip for the mean square displacement in the form 

<t) - — i' n^ o + 1- + -1^)1) • (A24) 

Here we have introduced the function , 

y{z) = / dy 

Jo y(i + r) 

= logz + 7E-^[e-"Ei(z) + e^Ei(-z)], (A25) 

where 7e = 0.577215665 is Euler's constant and Ei is the exponential integral [2^ . Using the expansion of the 
exponential integral for small argument, we obtain the expansion 



V{z) = -(logz + 7e) (cosh 2 - 1) - -[e-^ E ^ + E S^]' (^26) 

n—l n—1 

We see from this that for small z, 

^ -^^'(log ^ + 7E -^)-^z'(logz + 7E -§) + ••• . (A27) 

Note that V{Q) — 0, in agreement with the definition (jA25p . For large z, using the asymptotic formulas for the 
exponential integral, we obtain the asymptotic expansion, 

T/(z)^logz + 7E-^-^--^---- . (A28) 

With these results, we see that for very short times {t <C t) the mean square displacement again takes the form 
(|A13|) but now with the mean square velocity, given by 

^ [n^H - ^) + arccos ^ + 7^i»^ log ^ 

^ {--flogujQT + — ^ arccos-^}, (A29) 

Trm oji 2ujQ 

where second form is that for small relaxation time, ^ —t 1/t. In the Ohmic limit this is logarithmically divergent, 
so we have here a case where the single relaxation time makes a difference. On the other hand. 



7rma;i(ri2 — 7$! + uJq) 
arccos--^. (A30) 



Thus the Ohmic limit of (a;^) is finite. 
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d. Free particle 

The free particle corresponds to the absence of the oscillator potential, that is to the limit K ^ Q. For the Green 
function we can obtain this limit by setting = and wi = «7/2 in the expression (jAlSp . This gives 

In this free particle case the parameters J7 and 7 are given by the relations (|A15p with K = Q. These can then be 
inverted to give 



1 + y/l - AQr/m 1 - ^1 - ACr/m 
n= , 7= ^ . (A32) 

With this expression for the Green function, the high temperature form (|A2ip of the mean square displacement 
becomes 

2kT n3(i_g-y«)„^3(i_g-0«) 

'^'^ = jQin^ ^ 7^) ^ 

2kT , , 
= --(t ). (A33) 

C 7 

Note that at long time this increases linearly with time, consistent with ( IA12[) in the sense that (x^) — 00 for the 
free particle. On the other hand for short times, we get exactly the short time result (jAlSp with (i^) given by the 
equipartition form (|A22p . In other words, the short time behavior of the oscillator is that of the free particle. 
At zero temperature the result (jA24p becomes for the free particle. 



At very short times {t <C r) this takes the form (|A13p with now 



■^W = — 7^ 2 • ^-^^V 



(x^) = -^log^^-^log7r. (A35) 
7rm(si — 7) 7 TTTO 



At very long times (t 3> 7 we find 



2h 

^ —\ogCt/m. (A36) 



3. Weak coupling 

The coupling to the heat bath is measured by the function If this is small, the response function will be 

sharply peaked about — \J K/m^ the natural frequency of the oscillator. In the integral expression (|A5P for the 
Green function, one is therefore led to make the replacement /i(ti') —f i1{loo). The result is an expression of the Ohmic 
form (|MOl) with 

7 = -Re{/i(u;o)}. (A37) 

TO 

(The imaginary part gives a negligible contribution to wq-) Next, in the integral expression (jAlip for the mean square 
displacement, we make the same approximations with, in addition the replacement ?iwcoth^^ ftwo coth , to 
obtain 

s{t) ^ 2m (x^) I dt'G(t') 
Jo 

= 2{x^) {1- e-^*^^{cosLUit+ -^sinujit)}, (A38) 
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where 

>(.»>.|^eo.h^. (A39) 

This constitutes the weak couphng approximation (s^l • We should emphasize that this weak coupling approximation 
is not valid for the free particle, as should be clear from the above argument. Note, incidentally that for the Ohmic 
model this approximation is exact in the high temperature limit. The usual statement is that it is valid in the limit 
7 ^ but consideration of the exact results given above tells us that even in this limit the approximation is not 
correct for very short times {loqI <C 1) nor for very long times {-yt 3> 1). However, for all other times, the weak 
coupling approximation is very good for surprisingly large values of the coupling. As an illustration in Figure ?? we 
compare the mean square displacement at zero temperature as calculated first with the exact formula (|A24p and then 
with the weak coupling approximation. The parameters O/wq = 5 and 7/079 — 10/13 were chosen to exaggerate the 
difference. What we see is that the weak coupling approximation is surprisingly good, even for rather strong coupling. 

If one makes the further approximation of neglecting quantities of relative order j/ujo, one gets what at times is 
called the Weisskopf-Wigner approximation. In Eq. (|A38I| it would correspond to replacing ui —> ujo and dropping 
the second term after the exponential. This then would correspond exactly to what one obtains by solving the well 



known weak coupling master equation [31|. Because of this, in the literature the Weisskopf-Wigner approximation 
is often called the weak coupling approximation. The difference between the weak coupling approximations we have 
defined it and the Weisskopf-Wigner approximation is illustrated dramatically in Fig. ??. 

APPENDIX B: MATHEMATICAL FORMULAS 

Here we collect some formulas used in the evaluation of the various examples. These formulas are all simple and 
more or less well known. The first is the standard Gaussian integral, 

duexp{—^au^ + bu} = \—exp{^}. (Bl) 
2 \ a 2a 

The generalization to d dimensions takes the form 

1 (2n]'^^^ 1 
duexp(--u- A-u + B-u} = \ ' pxp{-R . A"! . B}. (B2) 

2 V dot A 2 

Here we have used dyadic notation, with A a positive definite symmetric matrix and B a vector (not necessarily real) 
in d dimensions. This generalization follows from the standard integral, using the fact that a symmetric matrix can 
be diagonalized by an orthogonal transformation. 

The second formula is the Baker-Campbell-Hausdorf formula. If A and B are a pair of operators (not necessarily 
Hermitian) whose commutator is a c-number, then 

gAgS ^ ^A+B^[A.B]/2 ^ ^B^A^[A,B] (gg) 

This formula is easily checked by expanding the exponentials in powers of their argument. A generalization of this 
theorem is the convenient formula 

e^g{B)=g{B+[A,B])e^, (B4) 

which holds for a general function g{B). Again, this can be verified by expanding g in powers of its argument. 

Finally, we have a couple of formulas based on the notion of a Gaussian variable. In general a set of operators 
(each with mean zero) is Gaussian if the expectation of a product of an odd number of the operators is zero while 
the product of an even number is equal to the sum of products of pair expectations, the sum being over all (2n — 1)!! 
pairings with the order within the pairs preserved. A Gaussian variable, e.g. x(i), is such a set with the members 
labeled with the time. Thus, for example, with an obvious shorthand, 

(1234) = (12) (34) + (13) (24) + (14) (23) . (B5) 

Note that within each pair the order is the same as the original order. A straightforward consequence of this Gaussian 
property is that for a Gaussian operator O, we have the formula 

(e''^)=e-Mc^>. (B6) 
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This is easily verifies by expanding the exponentials. Another convenient result is 



( 



1 



2^2 



2 



1 



a 



2 



\/27rcr2 



exp{ 



cxp{ 



2(a2 + (02)) 



}■ 



(B7) 



To obtain this result, form the Fourier transform with respect to the parameter a using the standard Gaussian integral 
(IBip and then form the expectation using (|B6p . 
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Figure 1. Probability distribution for a free particle in the Schro dinger cat state. P{x, 0) is the initial distribution. 
P{x,t) is the distribution at time t, while PQ{x,t) is the distribution obtained at time t by artificially setting the 
attenuation coefficient equal to zero. 

Figure 2. The motion of the wave packet center for the displaced ground state and for the coherent state, both for 
initial velocity zero. The displaced ground state motion is computed at zero temperature. The parameters chosen are 
7/1.^0 — 10/13 and Sl/o^o = 5 . 

Figure 3. The function A{t) for a coherent state pair. The parameters chosen are 7/wo = 10/13, fi/wo = 5. 

Figure 4. Constant density contours of the Wigner function for a squeezed state, shown in the plane of the 
dimensionless variables u = qj yj {x^) and v = p/m^y (i^)- The dashed circle corresponds to the equilibrium state, 
the state just before the initial squeeze as well as the state at long times. The contour marked (0) corresponds to the 
initial squeezed state. The contour marked (1/4) is that corresponding to a quarter period, while that marked (1/2) 
is that corresponding to a half period. The parameters chosen are 7/0^0 = 10/13, fJ/wo = 5. 

Figure 5. Second derivative of the Green function for the oscillator. Parameters for the single relaxation time model 
are 7/0^0 = 10/13 and fl/too = 5. 

Figure 6. Comparison of the exact and weak coupling expressions for the mean square displacement at zero 
temperature for the oscillator. The parameters chosen are j/lvq = 10/13 and fl/ujo = 5. 
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